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Kink defects in the 90° partial dislocation in silicon are studied using a linear-scaling density-matrix 
technique. The asymmetric core reconstruction plays a crucial role, generating at least four distinct 
kink species as well as soliton defects. The energies and migration barriers of these entities are 
calculated and compared with experiment. As a result of certain low-energy kinks, a peculiar 
alternation of the core reconstruction is predicted. We find the solitons to be remarkably mobile 
even at very low temperature, and propose that they mediate the kink relaxation dynamics. 
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The importance of dislocations in semiconductors 
hardly needs comment. In addition to being responsible 
for plastic behavior in general, dislocations occur com- 
monly at semiconductor interfaces where they can act 
as trapping and scattering centers for carriers. In sil- 
icon, the predominant slip system consists of 60°-edge 
and screw dislocations oriented along (110) and lying in 
a {111} slip plane. Both are known to dissociate into 
pairs of partial dislocations bounding a ribbon of stack- 
ing fault jl) . The resulting 90° or 30° partial dislocations 
are believed to have reconstructed cores, consistent with 
the low density of dangling bonds as observed by EPR 
measurements pj]. Since the dislocation motion occurs 
by nucleation and propagation of kinks along the dislo- 
cation line, a detailed understanding of the atomic-scale 
structure of the kinks is obviously of the greatest impor- 
tance. Unfortunately, experimental approaches have not 
proved capable of providing such an understanding. 

Until recently, the only theoretical methods capable 
of treating such problems were based on classical inter- 
atomic potentials. These are of questionable accuracy, 
and are generally unable to reproduce effects of intrinsic 
quantum-mechanical nature such as bond reconstruction 
and Peierls or Jahn- Teller symmetry breaking. For ex- 
ample, while the Stillinger- Weber || potential has been 
used to study the core reconstruction and kinks of the 
30° partial ||, it fails to reproduce the correct core re- 
construction of the 90° partial jjj. It is thus exciting to 
find that ab-initio methods are approaching the point of 
addressing some interesting questions about dislocations. 
Recent theoretical work has focused on such issues as the 
core reconstruction of the 90° ||,|| , and the elastic inter- 
action between dislocations of a dipolc in the shuffle || 
and glide M sets. One first-principles study has even 
been done on a kink barrier in the 30° partial, |?J but it 
assumed a kink structure previously proposed, and used 
a relatively small supercell. Cluster calculations on the 
90° partial have also been reported 01. However, a com- 



prehensive study of dislocation kink structure and dy- 
namics would require the use of very large supercells, for 
which the application of ab-initio techniques is still com- 
putationally prohibitive. Thus, there is a pressing need 
for the application of more efficient quantum-mechanics 
based methods to study the electronic and structural ex- 
citations in this system. 

In this Letter, we use a total-energy tight-binding 
(TBTE) description of the electronic and interatomic 
forces to carry out a detailed atomistic study of this 
kind for the 90° partial dislocation in silicon. The key to 
making the calculations tractable is our use of a "linear- 
scaling" or "O(N)" method of solution of the Schrodinger 
equation ||, enabling us to treat system sizes up to 
10 3 atoms easily on a workstation platform. We ver- 
ify that the dislocation core reconstructs with a sponta- 
neous symmetry breaking, and find that the "soliton" 
defect associated with the reversal of the core recon- 
struction is extremely mobile. We find that at least four 
distinct kink structures must be considered (labeled by 
the sense of the core reconstruction on either side) , and 
show that they can be classified as high- or low-energy 
kinks depending on whether or not they contain a dan- 
gling bond. Molecular-dynamics simulations as well as 
fully-relaxed static calculations are used to characterize 
formation energies, migration barriers, and kink-soliton 
reaction pathways. The picture that emerges is one in 
which the ground state is free of dangling bonds (even in 
the presence of kinks), and in which the solitons mediate 
the structural excitations and dynamics. 

We use the TBTE parameters of Kwon et al. (II]] , with 
a real-space cutoff of 6. 2 A on the range of the density 
matrix used in the O(N) method ||. We chose to work 
at a fixed electron chemical potential 0.4 eV above the 
valence-band edge, and thus all "energies" reported be- 
low are technically values of grand potential || . The nu- 
merical minimization of the 0(J\f) functional was carried 
out by the conjugate-gradient algorithm, with the inter- 



1 



FIG. 1. (a) Top view of the slip plane of a reconstructed 
90° partial dislocation. The shaded area indicates the stack- 
ing fault. Horizontal and vertical directions correspond to 
[110] and [112] directions, respectively, (b) Reconstruction 
defect, or soliton, where the core reconstruction changes ori- 
entation. 

nal line minimization performed exactly [ jTl| . Ground- 
state structures were computed by allowing all atomic 
coordinates to relax fully (forces less than 5 meV/A). 
The supercells used will be described below; all energies 
for defect (soliton and kink) structures are given with 
respect to a corresponding supercell containing defect- 
free (but reconstructed and fully relaxed) dislocations. 
Barrier energies were calculated by choosing a reaction 
coordinate and, for a series of values of this coordinate, 
computing the energy with this coordinate fixed and all 
others fully relaxed. Molecular dynamics runs were per- 
formed with a Verlet algorithm, the Nose thermostat [[L2| , 
and a time step of 2 fs. 

In Fig. 1(a), a top view of the atomic structure of the 
reconstructed 90° partial in its slip plane is shown. The 
shaded area indicates the stacking fault. The fourfold co- 
ordination of the atoms at the core is restored by bonding 
across the dislocation line. This reconstruction has been 
discussed by other authors ^||| and compared with an 
alternative "quasi-fivefold" bond reconstruction that pre- 
serves the mirror symmetries along the dislocation line. 
As a test of our TBTE-DM approach, we computed the 
energy difference between the two possibilities. In this 
calculation the supercell consisted of a 96-atom slab nor- 
mal to a (110) direction, containing two dislocations with 
opposite burgers vectors separated by a distance of 13.3 
A. We find the fourfold reconstruction to be 0.18 eV/A 
lower in energy than the fivefold one, in perfect agree- 
ment with the TB calculations in Ref. ||. Our results 
also compare favorably with the ab initio results in Ref. 
||, for which the energy difference between the two re- 
constructions is 0.23 eV/A. In our calculations, the re- 



is on the order of the energy associated with a dangling 
sp 3 orbital in silicon. 

An interesting question is whether a soliton can move 
easily along the dislocation. We computed an energy 
barrier of only 0.04 eV for the propagation of a soliton 
between two adjacent equilibrium positions, as indicated 
in Fig. 1(b). With such a small barrier, it might be ex- 
pected that the solitons would be extremely mobile even 
at very low temperatures. To test this, we performed 
a molecular dynamics simulation on a supercell having 
a soliton-antisoliton pair, initially separated by 9.6 A, 
on an otherwise defect-free partial dislocation. Remark- 
ably, at a temperature of only 50 K°, the solitons were 
indeed mobile and recombination of the pair took place 
after only 1.3 ps. Such highly mobile solitons play an 
interesting role in the relaxation of high-energy kinks, as 
explained below. 

A schematic view of the supercell we used for simulat- 
ing kinks and soliton-kink complexes is shown in Fig. 2. 
It contains a total of 864 atoms, corresponding to the 
96-atom slab repeated nine times along the dislocation 
line ([110] direction). The supercell vectors and the crys- 
talline directions are indicated. In Figs. 3(a)-(e) the lo- 
cal structures of the five different types of kinks are dis- 
played. The notation we chose to name each kink type is 
related to the orientation of the reconstructed bond, as 
one moves from left to right in each diagram in Fig. 3. 
For example, in Fig. 3(a) we denote the orientation at 
the left of the kink as "right" (R). Hence, we call this a 
right-left (RL) kink, the notation following accordingly 
for the other types. The above supercell was used to 
compute the energies for the LR and RL kinks. The lat- 
tice vector was staggered by twice the "kink vector" [[l4j , 
as shown in Fig. 2, in order to accommodate one RL and 
one LR kink in each dislocation. Fig. 3(f) shows a com- 
plex of a LR kink and a soliton. For the configurations in 
Figs. 3(c)-(f) the bonds on either side of the defect have 
the same orientation. For these, a supercell of 432 atoms 
was used (half of that shown in Fig. 2) , having one kink 
or complex in each dislocation, with a supercell vector 
staggered by one kink vector. 
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FIG. 2. Supercell used for kinks and kink-soliton com- 
plexes. Crystalline directions and supercell vectors indicated. 
Unit cell contains 864 atoms. 



Our results for the energies of the configurations in 
Fig. 3 are shown in Table I. It is seen that the LR and 
RL kinks are much lower in energy than all the others. 
(The LL and LL* kinks are found to be unstable against 
emission of a soliton, a mechanism that is discussed in 
detail in the next paragraph.) An inspection of Fig. 3 
shows that the LR and RL kinks arc fully reconstructed, 
with no dangling bonds. On the other hand, the high- 
energy RR kink and the unstable LL kinks all contain 
a dangling bond (note the three-fold coordinated atom 
at the core of each defect). This distinction is clearly 
responsible for most of the energy difference. On the 
basis of the energetics alone, we can thus conclude that 
any kinks which occur in the 90° partial dislocation in 
silicon will be almost entirely of type LR or RL. This 
means that the orientation of the core reconstruction is 
predicted to alternate from one inter-kink segment to the 
next along the dislocation. 

Moreover, we see that the RL + S and LR + S com- 
plexes (S = soliton) have lower energies than the RR and 
LL kinks. From the energetics in Table I we can write 
the following reaction equations: 



RR -> RL + S + 0.41 eV; 
RR — > LR + S + 0.36 eV; 
LL — ► LR + S (spontaneous); 
LL* — > RL + S (spontaneous). 



(1) 
(2) 
(3) 
(4) 



Reactions (1) and (2) above consists of the emission of a 
soliton by the RR kink, which turns into a LR or RL in 
the process. We estimated the energy barrier for these 
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FIG. 3. Core structure of various kinks and a kink-soliton 
complex. Kink notation relates to bond reconstruction on 
each side of the defect, and is explained in the text, (a) RL 
kink, (b) LR kink, (c) LL kink, (d) LL* kink, (e) RR kink, 
(f) LR kink + soliton. 
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TABLE I. Calculated formation energy for defects in the 
core of the 90° partial dislocation in silicon. Included also are 
the migration barriers for a soliton and the low energy kinks. 



Formation energy 
(eV) 



Migration barrier 
(eV) 



soliton 
LR kink 
RL kink 
LL kink a 
LL* kink a 
RR kink 

soliton + LR kink 
soliton + RL kink 



1.31 
0.50 
0.50 
1.74 
1.76 
2.04 
1.68 
1.63 



0.04 
1.87 
1.83 



"Approximate energy. Defect is unstable 

processes, obtaining the value of 0.05 eV. This result is 
not surprising in view of the low energy barrier for soliton 
motion. The same mechanism is involved in the instabil- 
ity of the LL and LL* kinks. Here the asymmetry of 
the local strain fields is enough to remove the barrier to 
emission of the soliton in one direction. 

We decided to put the above picture to test by per- 
forming a MD simulation on a 864-atom supercell con- 
taining four RR kinks, two in each dislocation. Thus it 
should be possible for the system to convert these kinks 
into alternating LR and RL kinks, as is to be expected 
from the energetics in Table I. At room temperature a 
RR kink in one of the dislocations emitted a soliton after 
0.2 ps, turning into a RL kink. This soliton propagated 
towards the other RR kink in the same dislocation and 
fused with it (converting it to a LR kink) after only 0.7 
ps. The latter process is equivalent to recombination of 
the propagating soliton with an anti-soliton "embedded" 
in the RR kink. 

In Table I, we also include the energy barriers to mo- 
tion of the RL and LR kinks. For metals, the formation 
of double kinks controls the rate of dislocation motion, 
and the energy barriers to kink propagation along the 
dislocation line are very small. The high values we ob- 
tained for silicon are a signature of the highly directional 
bonds in covalent semiconductors. In these systems, dis- 
location motion is believed to be controlled by the kink 
mobility. Recent experiments [|l5],[l6) have confirmed this 
picture, but some controversy still remains, and the role 
of impurities as obstacles to kink motion is yet to be fully 
explored. Here we concentrate on the important issue 
of the size of the kink formation energies and migration 
barriers, and make a connection with the experimental 
results. For the velocity of a gliding dislocation we have, 
from the theory of Ref. fTjJ 



v d oc exp[-(U k + W m )/kT] , 



(5) 



electron microscopy (TEM) or intermittent loading (IL) 
measurements [|l5|l6 range from 0.40 to 0.62 eV for Uk 
and 1.50 to 1.80 eV for W m . Our results from Table I, 
Uk = 0.50 eV and W m = 1.85 eV (average between LR 
and RL values), fall within the range of the experimental 
numbers. The theory in Ref. jl?]] can be used to cal- 
culate the Peierls stress of materials. When the above 
numbers are used for silicon, one obtains a value that is 
too low in comparison with results of high-stress mea- 
surements. Discrepancies are also found for quantities 
such as the velocity of steady-state motion of a disloca- 
tion under static load |l6fl . It has been argued in Ref. 
|]l6| that these discrepancies are to be assigned to point 
defects that would influence both the concentration of 
kinks and their migration barriers, by creating inhomo- 
geneities in the potential relief (the potential felt by a 
moving dislocation due to the periodicity of the lattice). 
The above comparison between experiment and our re- 
sults, which are valid for a homogeneous potential relief, 
reinforces the plausibility of this scenario. 

In summary, we investigated different structural and 
dynamical properties of the 90° partial dislocation in sil- 
icon. We verified that the core undergoes a reconstruc- 
tion that breaks the mirror symmetry along the dislo- 
cation line, as previously reported. This leads to the 
existence of solitons which are shown to be highly mobile 
along the dislocation core, and to a multiplicity of kinks 
whose stability is found to depend, in each case, on the 
reconstruction of dangling bonds at the core of the de- 
fect. We find that the high-energy kinks transform into 
low-energy ones by soliton emission. The energy barriers 
associated with this mechanism are small or absent in 
each case. The low-energy kinks are fully reconstructed 
with no dangling bonds, and impose an alternation of the 
orientation of the core reconstruction from one inter-kink 
segment to the next. The mobility of these low-energy 
kinks presumably determines the dislocation mobility in 
the glide plane. Our calculated formation energies and 
energy barriers for these kinks are in good agreement 
with available experimental estimates. 

During the final stages of preparation of this 
manuscript, we became aware of related work by Hansen 
et al. [p"8| , who consider some similar kink structures but 
with a much smaller interkink separation than was con- 
sidered here. 

This work was supported by NSF Grant DMR-91- 
15342. R. W. Nunes acknowledges the support from 
CNPq - Brazil. J. Bennetto acknowledges support of 
ONR Grant N00014-93-I-1097. 



where Uk is the formation energy of a kink and W m is 
the energy barrier for kink migration along the disloca- 
tion. The experimental estimates based on transmission 
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